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ABSTRACT 

We use published mid-IR and V-band flux ratios for images A and B of Q2237+0305 
obtained in September 1999 to demonstrate that the size of the mid-IR emission 
region has a scale comparable to or larger than the microlens Einstein Radius (ER) 
(~ 10 17 cm for microlensing by solar mass stars). Q2237+0305 has been monitored 
extensively in the R and V-bands for ^15 years. The variability record shows significant 
microlensing variability of the optical emission region, and has been used by several 
studies to demonstrate that the optical emission region is much smaller than the ER 
for solar-mass objects. For the majority of the monitoring history, the optical flux 
ratios have differed significantly from those predicted by macro-models. In contrast, 
recent observations in mid-IR show flux ratios similar to those measured in the radio, 
and to predictions of some lens models, implying that the mid-IR flux is emitted from 
a region that is at least 2 orders of magnitude larger than the optical emission region. 
We have calculated the likeli-hood of the observed mid-IR flux ratio as a function of 
mid-IR source size given the observed V-band flux ratio. The expected flux ratio for a 
source having dimensions of ~1 ER is a sensitive function of the macro model adopted. 
However we find that the probability of source size given the observed flux ratios is 
primarily sensitive to the ratio of the macro-model magnifications. Limits on the mid- 
IR source size can therefore be considered as a function of a one dimensional, rather 
than a 4 dimensional (two optical depths plus two shears) class of models. The majority 
of published macro models for Q2237+0305 yield a flux ratio for images B and A of 
0.8 - 1.1. By combining probabilities from the ratios A/B and C/D we infer that the 
diameter of a circular IR emission region is > 1ER with > 95% confidence. For sources 
of this size, other geometries, specifically an annular geometry, appropriate for a dusty 
torus, yield the same limit if the projected area rather than radius is considered. For 
microlensing by low-mass stars, this source size limit rules out non-thermal processes 
such as synchrotron as mechanisms for mid-IR emission. 

Key words: quasars - accretion - dust - individual: Q22347+0305 - gravitational 
lensing - microlensing 


1 INTRODUCTION 

Q2237+0305 (The Einstein Cross) was discovered in the CfA 
Redshift survey (Huchra et al. 1985). The object comprises a 
source quasar with a redshift of z = 1.695 that is gravitation¬ 
ally lensed by a foreground galaxy (z = 0.0394) producing 
4 images with separations of ~ 1". As a result of the prox¬ 
imity of the lensing galaxy, Q2237+0305 provides a unique 
opportunity to contrast dynamical measurements based on 
the geometry of the lensing with more traditional techniques 


(e.g. Rix, Schneider & Bahcall 1992, hereafter RSB92). In 
addition, the close proximity results in a short time-scale 
for microlensing, and a large projected microlensing length- 
scale (with respect to the source). Each of the 4 images are 
observed through the bulge of a galaxy which has an opti¬ 
cal depth in stars that is of order unity (e.g. Kent & Falco 
1988, hereafter KF88; Schneider et al. 1988, hereafter S88; 
Schmidt, Webster & Lewis 1998, hereafter SWL98). This 
results in a high probability for microlensing. 

Microlensing is indicated either by independent tem- 


© 1999 RAS 


2 J. S. B. Wyithe et al. 


poral variability of image fluxes, or by variation in colour 
between images at a single epoch (separated by the macro¬ 
image delay). Large, and rapid variation of the continuum 
flux is found in the variability record for Q2237+0305 (Irwin 
et al. 1989; Corrigan et al. 1991; 0stensen et al. 1996; Woz- 
niak et al. 2000a,b). This variation has been used to argue 
that the optical emission region must be significantly smaller 
than the microlens Einstein Radius (ER), and therefore 
the typical scale of the caustic structure (e.g. Wambsganss, 
Paczynski & Schneider 1990; Wyithe, Webster, Turner & 
Mortlock 2000; Yonehara 2001). During a caustic crossing, 
a small source exhibits colour variability if the emission 
spectrum is scale dependent (e.g. Wambsganss & Paczynski 
1991; Fluke & Webster 1999). Evidence of this effect from 
broad band observations of a microlensing event was pre¬ 
sented by Corrigan et al. (1991). Furthermore, if the quasar 
emits in one wave-band at a scale much smaller than the 
ER, and on a scale larger than the ER in another band, 
then colour change may be seen in two random observations 
of a single image (particularly if the observations straddle a 
caustic crossing event), or between two different images, as a 
result of magnification of the smaller source. Ground based 
observations have confirmed differential amplification of the 
emission region. Lewis et al. (1998) determined the ratios of 
emission line equivalent widths relative to one image. They 
show ( i ) the ratios remain fairly constant for one image from 
line to line, suggesting that the sizes of the emission regions 
for the lines are not greatly different, ( ii ) that the ratios 
vary from image to image for a single epoch by a factor of 
~ 2.5, and (in') that the ratio for a single image varies as 
a function of time, i.e. as a result of a microlensing event. 
These results are consistent with earlier results of Fillipenko 
(1989) who measured a ~ 25% difference in the width of the 
Mgll lines between the A and B images. 

The CIII] line, produced by extended broad line regions 
has been measured in an attempt to find the flux ratios ( R ) 
using emission scales beyond the influence of microlensing 
(Yee & De Robertis 1992; Racine 1992; Fitte & Adam 1994; 
Saust 1994; Lewis et al. 1998). Mediavilla et al. (1998) ob¬ 
served the CIII] line in Q2237+0305 using two-dimensional 
spectroscopy, and found an arc (image of the extended nar¬ 
row line region of the source) extending around three of the 
images, indicating a very extended region of emission. On 
the other hand, measurement of R using the CIII] line is sub¬ 
ject to differential extinction, and uncertainties in contin¬ 
uum subtraction. To avoid the effects of extinction, and pos¬ 
sibly microlensing Falco et al. (1996) imaged Q2237+0305 
at 3.6cm, finding flux ratios similar to those inferred from 
CIII]. In addition, Q2237+0305 has also been imaged in the 
Ultraviolet (Blanton, Turner & Wambsganss 1998), and in 
X-rays (Wambsganss, Brunner, Schindler & Falco 1999). 

Many models have been proposed for the projected lens 
mass distribution based on observations of the lensed images 
of Q2237+0305 (e.g. KF 88 ; S 88 ; Kochanek 1991, hereafter 
K91; Wambsganss & Paczynski 1994, hereafter WP94; Witt, 
Mao & Schechter 1995, hereafter WMS95; SWL98; Chae, 
Turnshek & Khersonsky 1998, hereafter CTK98). The ma¬ 
jority of these predict a flux ratio for images B and A of 
Rba ~ 0.80 — 1 . 1 , consistent with the ratio measured in the 
radio of Rba = 1.1 ± 0.3 (Falco et al. 1996). However, over 
the monitoring history, the optical light-curve shows vari¬ 
ations in Rba between ~ 0.2 and ~ 1.0. The discrepancy 


Table 1. Values of the optical depth (re) and the magnitude of 
the shear ( It I )f° r each simulation. The number of stars (IV*), the 
theoretical mean magnification ((jiih) = |(1 — re) 2 — 7 2 | — 1 ), and 
the mean magnification of the map f {//}) are given in each case. 


model 

K, 

M 

TV* 

(iHh) 

(v) 

i 

0.30 

0.35 

8374 

2.72 

2.75 

ii 

0.30 

0.40 

11022 

3.03 

3.05 

Hi 

0.30 

0.45 

15471 

3.48 

3.50 

iv 

0.30 

0.50 

23724 

4.17 

4.26 

V 

0.35 

0.35 

13253 

3.33 

3.48 

vi 

0.35 

0.40 

18501 

3.80 

3.95 

vii 

0.35 

0.45 

28262 

4.55 

4.79 

viii 

0.40 

0.35 

21709 

4.21 

4.29 

ix 

0.40 

0.40 

32993 

5.00 

4.99 

X 

0.45 

0.35 

37964 

5.55 

5.61 

xi 

0.69 

0.71 

16915 

2.45 

2.40 


is attributed primarily to microlensing, and firmly demon¬ 
strates that the optical flux ratios cannot be used as model 
constraints (e.g. S 88 ; KF 88 ; K91; WP94). Agol, Jones & 
Blaes (2000) (hereafter AJB00) have found a mid-IR B:A 
flux ratio of Rba ~ 1.1. This ratio is consistent with obser¬ 
vations in the radio (Falco et al. 1996). AJB00 interpreted 
their results as evidence for an extended region of mid-IR 
emission, with dimensions larger than the microlens Einstein 
Radius. 

In this paper we use microlensing models to calculate 
distributions of flux ratios for sources with different sizes and 
intensity profiles, and hence derive quantitative limits on the 
scale of the mid-IR emission. In Secs. ^ and we describe 
the microlensing models, and summarise published macro¬ 
models for Q2237+0305. Sec. [jj discusses the methods used 
to infer the mid-IR source size from the observed optical and 
mid-IR flux ratios, and the source size limits implied by the 
published macromodels. Initially we restrict our discussion 
to the particular case of the flux ratio between images B and 
A. However we present results based on all image ratios in 
Sec. 4.3. In the conclusion we mention some implications for 
quasar physics. 


2 MICROLENSING MODELS 


Throughout the paper, standard notation for gravitational 
lensing is used. The Einstein Radius is defined to be the 
radius inside which the mean surface mass density is equal 
to the critical density. The Einstein Radius r/o of a 1 M@ star 
projected into the source plane is 


Vo 


DdsD a 4GMq 


D d 


(1) 


where G is the gravitational constant, c is the speed of light, 
and D d , D s and D ds are the angular diameter distances of 
the lens, source and from the lens to the source respectively. 
The normalised shear is denoted by 7 , and the convergence 
or optical depth by re. The model for gravitational microlens¬ 
ing consists of a very large sheet of point masses that simu¬ 
lates the section of the galaxy along the image line-of-sight, 
together with a shear term that includes the perturbing ef¬ 
fect of the mass distribution of the lensing galaxy as a whole. 
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The normalised lens equation for a field of point masses with 
an applied shear in terms of these quantities is 


y = 


1-7 ° 

0 1+7 J Z —— a?| 2 

' 7 = 1 


( 2 ) 


Here x and y are the normalised image and source positions 
respectively, and the x J t and m J are the normalised positions 
and masses of the microlenses. We have assumed that the 
entire optical depth is in compact objects, and that all mi¬ 
crolenses have a common mass. The ray-tracing method (e.g. 
Kayser, Refsdal & Stabcll 1986; Schneider & Weiss 1987; 
Wambsganss, Paczynski & Katz 1990) was used to com¬ 
pute the magnification maps. Each map produced covered 
an area of 100 ry o x lOOr/o, and was computed at a resolution 
of 2000 x 2000 pixels. The maps therefore describe the mag¬ 
nification distributions for a dynamic range in source size of 
~ 500. 100 rays were traced per unlensed pixel. The region 
of the lens plane in which image solutions need to be found 
to ensure that 99% of the total macro-image flux is recov¬ 
ered from a source point was described by Katz, Balbus & 
Paczynski (1986). The union of areas in the lens plane cor¬ 
responding to the flux collection area of each point in the 
area of the magnification is known as the shooting region, 
the method for determining the dimensions of which is de¬ 
scribed in Lewis & Irwin (1995), and Wyithe & Webster 
(1999). The radius of the disc of point masses was chosen to 
be 1.2 times that required to cover this shooting region. 

Tab. [I] shows the parameters k and 7 for each of the 
11 magnification maps produced, along with the number of 
stars used in the model, the theoretical mean magnification, 
and the mean magnification of the map. Models i — x cover 
the k — 7 parameter space near the values estimated for im¬ 
ages A and B by SWL98. A model (corresponding to image 
C of SWL98) with larger k and 7 (model xi ) is included 
for comparison, and as an indicator of the sensitivity of our 
results to the macro-parameters. 


3 MACRO-MODELS FOR Q2237+0305 

The quadruple quasar Q2237+0305 has been the target of 
numerous studies. The authors have attempted to explain 
the geometry of the system via gravitational lensing using 
various models. S 88 and later RSB92 assumed that mass 
in the lensing galaxy was proportional to observed light. 
Other approaches parameterise the lens using various sim¬ 
ple potentials having either intrinsic or external quadrupoles 
(KF 88 ; K91; WP94; WMS95; SWL98; CTK98). Constraints 
on the lens model come from the position of the galaxy and 
of the 4 quasar images. In addition, the flux ratios poten¬ 
tially provide 3 additional constraints. However, these are 
difficult to measure from optical images due principally to 
the effects of differential absorption, and microlensing. KF 88 
used the available optical flux ratios, while noting that the 
mismatch of the flux for images C and D is the dominant 
contribution to the error in their fit. Optical flux ratios were 
also used by S 88 and RSB92, though the minimisation quan¬ 
tity was weighted 9:1 in favour of image positions. More 
reliable flux ratios, from the CIII] line (Racine 1992) were 
employed by WMS95, and more recently, CTK98 have incor¬ 
porated the 3.6cm fluxes measured by Falco et al. (1996). 
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Figure 2. The cumulative probability of the mid-IR flux ratio 
= /<b Il l A as a function of mid-IR source size (7/ /+) for a 
variety of Gaussian source sizes Sip = 2cr = 0.4, 1.6. 0.4. 25. fir/,, 
(left to right). The parameters for both images A and B were 
K = 0.35, 7 = 0.40. Plots are shown for a uniform optical 
source: Sopt — 0.025 r/o (single pixel, solid lines), and 2 Gaussian 
sources: Sopt = 2 a = (I.O 077 (dashed lines) and Sopt = 2a = 
0. 1 r] t , (dot-dashed lines). The thick line represents the prior for 
optical flux ratio. 


Other studies fitted only image positions. All these stud¬ 
ies have reproduced the best measured positions (currently 
those of Blanton et al. 1998) of the time to within (or near 
to) the quoted error. However flux ratios, particularly in the 
optical have been less successful. A major conclusion of the 
parametric studies is that the image configuration (and flux 
ratios) can be reproduced by families of models, and that 
the image geometry alone is insufficient to uniquely com¬ 
pute the optical depth and shear at the location of each 
image. Hence parametric lens models do not uniquely define 
the individual image magnifications nor the total magnifi¬ 
cation. Indeed, different potentials predict total magnifica¬ 
tions ranging from a few to a few hundred (WP94; WMS95; 
CTK98). On the other hand, the total mass inside the Ein¬ 
stein ring is consistently measured at ~ 1.5 x 10 10 Mq/i^ 5 1 
(RSB92; WP94; CTK98). In addition, the large variation 
in total magnification is not found in the magnification ra¬ 
tios. Most published models have flux ratios for images B 
and A of ~ 0.80 — 1.1, and the degenerate model families of 
WP94, WMS95 and CTK98 all have Rba ~ 0.85 — 0.90. In 
a more physically motivated parametric model, SWL98 in¬ 
cluded a bar component in the lensing galaxy to correct for 
the apparent misalignment between the model major axis 
and observed galactic orientation. SWL98 broke the degen¬ 
eracy in their family of models (analogous to the shear - 
power-law degeneracy discussed by WP94) through the use 
of an elliptical mass model having the observed eccentricity 
of the light distribution. The collection of models from the 
literature described in this section are summarised in Tab. ^ 
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Table 2. Summary of published macro-models from the following authors Schneider et al. (1988) (S88) Kent & Falco (1988) (KF88), 
Kochanek (1991) (K88), Rix, Schneider & Bahcall (1992) (RSB92), Wambsganss & Paczynski (1994) (WP94), Witt, Mao & Schechter 
(1995) (WMS95), Schmidt, Webster & Lewis (1998) (SWL98), Chae, Turnshek & Khersonsky (1998) (CTK98). The labels corresponding 
to Figs. |)| [(J and |jj are shown in column ii. Where a paper presents multiple models, they are labelled in column iii using the authors 
nomenclature. The table presents values for the optical depth and shear (where available), the theoretical magnification for each image, 
and the magnification ratios. 


Author 

label 

model 


A 

171 

ImI 

K 

B 

171 

M 


C 

M 

M 


D 

M 

ImI 

r BA 

R 

r AC 

r DA 

S88 

S 


0.36 

0.44 

4.6 

0.45 

0.28 

4.5 

0.88 

0.53 

3.8 

0.61 

0.66 

3.6 

0.98 

1.20 

0.78 

KF88 

KF 

De Vauc. 



3.34 



2.64 



1.14 



2.11 

0.88 

2.63 

0.63 

K91 

K1 

1 



10.0 



6.8 



7.20 



8.10 

0.68 

1.43 

0.81 


K2 

2 



4.3 



3.83 



2.11 



5.07 

0.89 

2.04 

1.18 


K3 

3 



8.3 



7.30 



4.65 



8.88 

0.88 

1.78 

1.07 


K4 

4 



2.3 



2.05 



1.08 



2.23 

0.89 

2.13 

0.97 


K5 

5 



2.6 



1.72 



1.66 



1.72 

0.66 

1.56 

0.66 

RSB92 

RSB1 

1 

0.42 

0.42 

5.9 

0.41 

0.41 

6.1 

0.71 

0.63 

3.1 

0.65 

0.60 

4.1 

1.03 

1.92 

0.69 


RSB2 

2 

0.39 

0.46 

5.8 

0.38 

0.46 

6.3 

0.65 

0.67 

3.0 

0.59 

0.67 

3.6 

1.09 

1.96 

0.62 


RSB2a 

2a 

0.41 

0.47 

6.2 

0.38 

0.43 

5.1 

0.65 

0.68 

2.8 

0.59 

0.56 

5.1 

0.82 

2.22 

0.82 


RSB3 

3 

0.42 

0.41 

4.5 

0.37 

0.48 

4.4 

0.68 

0.66 

3.1 

0.65 

0.67 

3.6 

0.98 

1.45 

0.80 

WP94 

WP1 

0 = 0.0 

0.00 

0.77 

2.43 

0.00 

0.73 

2.12 

0.00 

1.36 

1.16 

0.00 

1.19 

2.38 

0.88 

2.08 

0.98 


WP2 

0 = 1.0 

0.47 

0.41 

8.98 

0.46 

0.40 

7.71 

0.56 

0.63 

5.04 

0.51 

0.58 

9.64 

0.86 

1.79 

1.07 


WP3 

0 = 1.7 

0.84 

0.13 

103.6 

0.84 

0.13 

87.80 

0.88 

0.17 

66.0 

0.86 

0.17 

118.5 

0.85 

1.56 

1.14 

WMS95 

WMS1 

0 = 0.5 

0.23 

0.60 

4.3 

0.22 

0.60 

3.8 

0.29 

0.96 

2.4 

0.26 

0.88 

4.4 

0.88 

1.79 

1.03 


WMS2 

0 = 1.0 

0.47 

0.42 

9.4 

0.46 

0.41 

8.1 

0.55 

0.62 

5.6 

0.52 

0.58 

10.0 

0.86 

1.56 

1.16 


WMS3 

0 = 1.5 

0.73 

0.22 

36.5 

0.72 

0.21 

31.3 

0.79 

0.30 

23.3 

0.76 

0.29 

40.7 

0.86 

1.56 

1.16 

CTK98 

CTK1 

e = 0.3 



3.69 



3.32 



1.86 



3.14 

0.9 

2.00 

0.85 


CTK2 

e = 0.2 



9.03 



8.15 



5.42 



9.03 

0.9 

1.67 

1.00 


CTK3 

e = 0.1 



27.4 



24.66 



17.81 



30.14 

0.9 

1.54 

1.10 

SWL98 

SWL1 

A = 0.5 

0.36 

0.40 

3.90 

0.36 

0.42 

4.39 

0.69 

0.71 

2.42 

0.59 

0.61 

4.91 

1.12 

1.63 

1.26 


SWL2 

A = 1.0 



4.05 



4.37 



2.47 



5.10 

1.08 

1.64 

1.26 


SWL3 

A = 2.0 



4.18 



4.34 



4.35 



2.51 

1.04 

1.67 

1.24 


Table 3. Values for the fractions of total flux in each image, and the flux ratios in the mid-IR (AJBOO) and V-bands (Wozniak et al. 
2000a,2000b). The V-band fluxes have been de-reddened (AJBOO), and the dominant contribution to the quoted error comes from this 
procedure. 





Fraction of Flux 


Wavelength 

Date (1999) 

A 

B 

C 

D 

8.9 & 11.7/im 

28th July Sz 21st September 

0.27T.02 

0.30T.02 

0.16T.02 

0.27T.02 

V-band 

1st August 

0.39T.08 

0 .11T.02 

0.41T.14 

0.10T.03 

V-band 

26th September 

0.46T.10 

0 .12T.02 

0.32T.11 

0.10T.03 


Wavelength 

Date (1999) 

nobs 

k ba 

nobs 

k ac 

Flux Ratio 

nobs nobs 

^DA ^BC 

nobs 
^ DB 

nobs 
^DC 

8.9 & 11.7/rm 
V-band 
V-band 

28th July & 21st September 
1st August 

26th September 

l.ii.i 

0.28T.08 

0.26T.07 

1.7T.3 

0.95T.37 

1.43T.55 

1 .0T.1 

0.26T.09 

0.22T.08 

1.9T.3 

0.27T.09 

0.38T.15 

0.9T.1 

0.91T.31 

0.83T.29 

1.7T.2 

0.24T.10 

0.31T.13 


4 MID-IR SOURCE SIZE LIMITS FROM FLUX 
RATIOS 

The magnification of an extended microlensed source is 
calculated by convolving the source intensity profile with 
a magnification map. Fig. |l] shows a magnification map 
(k = 0.35, 7 = 0.40) convolved with discs having radii of 
S'=0.15, 0.60, 2.10 and 7.35rj 0 . Circles corresponding to these 
sizes are plotted in the top-right corner of each map. The 
grey-scale bar describes the magnification. Fig. [i] demon¬ 
strates that even a source several ER in extent can be mi¬ 
crolensed by an appreciable factor (see for example Refsdal 
& Stabell 1993). In addition, there is a correlation between 
the caustic regions of highest density and the higher mag¬ 
nification regions for the larger sources. However the figures 
demonstrate that there is significant loss of correlation be¬ 
tween the magnifications of a source which is <C r/ 0 and 
of larger co-positional sources. In particular, larger sources 
show light-curves with lower event amplitudes and longer 


event durations (see the light-curves and correlation func¬ 
tions for various source sizes in Wambsganss, Paczynski & 
Katz 1990). 

The cumulative probability for the mid-IR flux ratio 
Rba — i f° r different mid-IR source sizes Sir given the 
measured optical ratio is 

/> OO 

r> / rRK . f~, r. \ / i t jO PT,obs / t jO PT,obs\ 

Pr{Rba< - \Sir,Sopt)= / dR BA P(R ba ) 

Jo 

xP r (R i b r a < —\Sir, Sopt, R BA T ’ obs )• (3) 

PA 

P(R < BA T '° bS ) i s the probability for the observed optical flux 
ratios, which we have assumed to be Gaussian with a mean 
and halfwidth (a) equal to the observed value and uncer¬ 
tainty. The integral over R BA T ’ ob3 was performed via Monte- 
Carlo. The observed optical flux fractions and ratios are 
summarised in Tab. |j| In the V-band (Wozniak et al. 2000b), 
the de-reddened flux ratios (AJBOO) were R < ^ A T,obs = 
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Figure 1. Convolution of a magnification map with source discs of 4 different sizes (shown by the circles in the upper right corner of 
each panel, the radii are 5=0.15, 0.60, 2.10 and 7.35r) 0 ). The grey-scale bar shows the corresponding magnification. The macrolensing 
parameters were k = 0.35, 7 = 0.40. The axes are in units of r) 0 - 

0.28 =t 0.08 on the 1st August 1999 and 0.26 =t 0.07 on the 
26th September 1999. In the remainder of this paper we 
use the V-band fluxes from the 26th of September 1999. 

Fig. 1 shows the cumulative probability for mid-IR flux ratio 
Pr(Rba < R\Sir, Sopt) for four different mid-IR source 
sizes Sir given three different optical source sizes Sopt in 


each case. The mid-IR sources were assumed to be Gaussian 
with half width a: Sir = 2cr = 0.4,1.6, 6.4, 25.6?7 0 (left to 
right). The optical sources were uniform: Sopt = 0.025?7o 
(single pixel), and Gaussian: Sopt = 2 a = 0.05»7 o (dashed 
lines) and Sopt = 2 a — O.lr/o (dot-dashed lines). The pa¬ 
rameters for both images A and B were k = 0.35, 7 = 0.40. 
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Figure 3. Left: The cumulative probability for the IR source size P(Sir < S\Sopt )? and Right: The differential probability for Sir 
given a uniform optical source: Sopt — 0.02577 o (single pixel, solid lines), and 2 Gaussian sources: Sopt = 2a = 0.05r) o (dashed lines) 
and Sopt = 2a = 0.1r] o (dot-dashed lines). The parameters for both images A and B were k = 0.35, 7 = 0.40. The mid-IR source was 
also assumed to be Gaussian with halfwidth a = Sir/2. 






Figure 4. The probability that the IR source size is greater than lr]o (1 — P (Sir < 1t]o\Sopt)) as a function of mean magnification. 
The parameters for both images A and B were the same in each case, the macro parameters were from models i — x (dark dots) and 
model xi (open circle). The optical source was uniform with Sopt — 0.025t7 o (single pixel). The mid-IR source was also assumed to be 
Gaussian with halfwidth a = Sir/2. 
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Figure 5. Plot of the probability that Sir > 1.0 r/ 0 as a function of mean flux ratio, Sopt = 0.025r; o (single pixel). The larger markers 
indicate models that include the large k large 7 for one of the images (model xi). The labels correspond to published macro-models (see 
Tab. ^j) and are plotted at the appropriate Rab- 


The thick line shows p(R < j^ I f r ’ oba ). Small mid-IR sources have 
distributions of R 2 Ja with means close to the observed op¬ 
tical ratio, while larger sources have distributions which are 
narrow, and symmetric about the macro-model flux ratio 
Rba (1 for the example in Fig. ^|). The symmetry signifies a 
loss of correlation between the flux ratio of the large mid-IR 
source and the smaller V-band source. Mid-IR source sizes 
of ~ 1 r) 0 show intermediate behaviour. R T ba for a given Sir 
is typically larger for smaller assumed Sopt due to the lower 
correlation of mid-IR and optical fluxes. 

Source size limits are imposed by the difference in the 
flux ratios for the optical source (which is known to be -C t] 0 ) 
and the mid-IR source. Given the observed V-band flux ratio 
(R < BA T ’° ba ) 1 the probability that the mid-IR flux ratio {R’ba) 
of two images B and A is smaller than the observed value 
R I BA° bS i s n °t a monotonic function of the mid-IR source 
size Sir for R BA greater than Rba (the macro-model flux 
ratio). However, the R ! ba f° r sources with Sir ~ r/ 0 are 
quite uncorrelated with R'gfT■ 

For each mid-IR source size we construct the likelihood 

7 jdir (Rba\Sir> Sopt) = 

dR BA 

at I I j-fIR.obs dPfl , J-.IR j-,IR,obs\ / r-,IR,obs\ / ,\ 

N J dR B A -Jfti Rba = Rba )p(R B A )• ( 4 ) 

P(R I BA ° ba ) the probability for the observed mid-IR flux 
ratio, which we have assumed to be Gaussian with a mean 
and halfwidtli (<r) equal to the observed value and uncer¬ 
tainty, and is from Eqn. The observed flux frac¬ 

tions and ratios are summarised in Tab. N. In the mid-IR 
(both 8.9 and 11.7 pm), AJB00 quote R B A° ba = 1.1 ± 0.1 
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for observations close to the V-band observations of Woz- 
niak et al. (2000) (28th July and 21st September). The in¬ 
finity in the integral over R was approximated by a 
value > 2R I J^ obs . Eqn. ^ is a function of the optical source 
size (Sopt) assumed. Sopt is thought to be ^ 10” 2 ri o 
from a variety of arguments (e.g. Wambsganss, Paczynski & 
Schneider 1990; Wyithe, Webster, Turner & Mortlock 2000). 
Rather than convolve with a further probability for Sopt, 
we have computed the likelihood (Eqn. [lj) for several values 
of Sopt < 0.1r; o . 

To construct probabilities for source size we have com¬ 
bined Eqn. with uniform logarithmic Bayesian priors 

where Sir < Smax 
otherwise (5) 

for the unknown mid-IR source size Sir: 


dlog(SiR) 


oc 1 


= 0 


Ps(S < Sir\Sopt) = 
r s ia 


N 


dS‘ 


/ dPp 1 


dPih 


dS 


dR T B R A 


( Rba\S Sopt). 


( 6 ) 


N is a normalising constant. Fig. shows the cumulative dis¬ 
tribution Ps(Sir < S\Sopt) as well as the differential prob¬ 
ability for source size corresponding to the example in Fig. ^]. 
The distributions are plotted assuming 4 different upper cut¬ 
offs in the prior for source size, Smax = 20,100, 500 and 

2500 r/ 0 . Smax is known to be smaller than ~ 25000 
since it is not resolved in the imaging data (AJB00). For 
this macro model, the source is larger than Sir = I77 0 with 
greater than 97% ( Sopt = 0.025r7 o ) confidence. The prob- 
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Figure 6. Plots of source size limits as a function of mean flux ratio. Results for three assumptions of optical size are shown; A 
uniform source: Sopt = 0.025r; o (single pixel, dots), and two Gaussian sources: Sopt = 2 <r = O.Oor/o (crosses) and Sopt = 2 cr = 0.1 rj 0 
(squares). Panels (a)-(c): Lower limits (Sir >0.25, 0.5 and 1.0?jo). Panels (d)-(f): Upper limits ( Sir <1.0, 5.0 and lO.Orjo). The larger 
markers indicate models that include the large k large 7 for one of the images (model xi). Note that panel (d) = one minus panel (c). 


ability shows minimal dependence on Sopt- In this case 
the measured R I ^ >bs differs from Rba only by the ob¬ 
servational uncertainty, and therefore carries no informa¬ 
tion on how large the source can be. As a result the up¬ 
per limits are entirely dependent on the upper cutoff for 
the prior assumed. The upper cutoff is the parameter to 
which the probabilities are most sensitive. Cumulative dis¬ 
tributions Ps(Sir < S\Sopt) were computed where images 


A and B are assumed to have identical parameters, hence 
Rba = 1 in each case. The confidence that Sir > 1 r / 0 is 
plotted as a function of the macro-model magnification in 
Fig. 0 assuming 4 different cutoffs Smax, the open circle de¬ 
notes model xi. There is some spread in these values, (e.g. 
Ps ~ 0.94 — 0.97 for Smax = 20 r/ 0 and Ps ~ 0.98 — 0.99 for 
Smax = 500t? o ), and higher mean model magnification tends 
to result in slightly tighter constraints. 
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Ps{S < Sir\Sopt ) was computed for all 121 combi¬ 
nations of the macro-parameters i — xi. Fig. shows the 
confidence that Sir > lr/o as a function of the macro-model 
flux ratio. Probabilities corresponding to Smax = 100 r/o are 
plotted (similar results are obtained for other values of Smax- 
Different assumptions for macro-models provide different 
limits since large sources have a flux ratio that is corre¬ 
lated more with the macro-model flux-ratio, and less with 
the optical flux ratio. Fig. |E| demonstrates a correlation be¬ 
tween probability and Rba that has surprisingly little scat¬ 
ter. When considering source size limits determined from 
flux ratios, we can therefore approximately parameterise dif¬ 
ferent macro models by their predicted flux ratio. The larger 
markers in Fig. denote those probabilities computed as¬ 
suming one or both image macro-models to be model xi. 
These lie below the other points, illustrating the generalisa¬ 
tion of the correlation in Fig. |f|. However, because the limits 
obtained from macro-models that include model xi follow 
nearly the same correlation as combinations of models i — x, 
the correlation may be valid for combinations of ka/7a and 
/ys other than those investigated. 

Distributions of flux ratios vary with macro-model for 
fixed Rba- Therefore the correlation in Fig. ^ is presumably 
due to the imposition of the optical flux ratio. The following 
approximation demonstrates why this should be the case. 
Consider two microlensed images b and a. Suppose that the 
magnification map for image a is uniform, while the mag¬ 
nification map for image b contains a single caustic but is 
otherwise uniform. Let the magnification of image a be 
Furthermore, let the magnification of image b be comprised 
of the magnification due to critical images associated with 
the caustic fi b in addition to the magnification from non- 
critical images fi b . The resulting optical magnification ratio 
is 


OPT,o , OPT,c OPT,c 

nOPT _ Mb + Mb p Mb_ 

ctba OPT,o ~ rLba ' OPT,o > 

f-^a 


(7) 


where Rb a is the theoretical flux ratio. Similarly, the larger 
mid-IR source size has a magnification ratio 


IR,c 

R T? i ftb 
ttba ~ Mba H-7^7 • 

/^■a 


( 8 ) 


Now /.ta PT '° = nl R '° = Ma by construction, so Rif —> R ba 
for large mid-IR sources. Consider a different model having 
both an increased fia and n b PT ’ 0 while keeping R ba con¬ 
stant. To maintain the observed R b cf T , lxf PT,c must be in¬ 
creased, either by increasing the caustic strength (flux fac¬ 
tor, see Witt (1990)), or by moving the optical source closer 
to the caustic. In the former case n b R ' C is increased by the 
same fraction as /r® PT ’ c , thus the same source size is re¬ 
quired to produce the observed flux ratio for different models 
having the same R ba . In the second case, shifting the larger 
mid-IR source relative to the caustic has a smaller effect on 
fj, b R ’ c than on [if FT ’ c , thus increasing [if PT ’° and /i® PT, ° 
slightly lowers the resulting R b f. A larger mid-IR source 
is therefore required to reproduce the observed mid-IR flux 
ratio, and hence the confidence for each (upper) limit is in¬ 
creased. The combination of these two effects suggests an 
explanation for both the trend (of tighter constraints with 
higher mean model magnification) shown in Fig. [i], and the 
correlation seen in Fig. H. 

The left-hand panels of Fig. |b| show the probabili- 
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ties (lower limits) that Sir > 0.257/ 0 , Sir > 0.5r/ o and 
Sir > 1.0 r/o as a function of the macro-model flux ratio 
Rba- Similarly, the right-hand 3 panels of Fig. |g| show the 
confidences (upper limits) that Sir < 1.0r/ o , Sir < 5.0 r] 0 
and Sir < 10.0?7o. Points are shown representing the three 
assumptions for the optical source size, a uniform source 
Sopt = 0.025 7] 0 (one pixel, dots), and 2 Gaussian sources 
Sopt = 2 a = 0.05r7 o (crosses) and Sopt = 2a = 0.1 r/ 0 
(squares). The confidences show some dependence on Sopt- 
Stronger limits are obtained if the optical source is larger 
since it is then not subject to such large fluctuations. In 
the remainder of this paper, we quote limits based on the 
assumption of Sopt = 0.025 r] 0 . 


4.1 IR-source size limits 

The statistics of flux ratios, particularly for large sources are 
sensitive to the macro-model assumed. Rather than select a 
particular model, or attempt to assign relative weights to 
limits found using different published models, we have com¬ 
puted the statistics for an ensemble of microlensing models 
for the images A and B (models i — xi in Tab. 111). These 
models cover the k — 7 parameter space near the values es¬ 
timated for images A and B by SWL98 (models i — x), as 
well as one with larger k and 7 (corresponding to image C of 
SWL98, model xi). The models have mean magnifications 
ranging from 2.72-5.56. Combinations of these models for 
images A and B of Q2237+0305 yield macro-model flux ra¬ 
tios Rba that differ by a factor of > 2 in both directions 
from the value observed in the mid-IR by AJB00. Most of 
the published macro-models (Tab. |i|) lie within the covered 
region, both in terms of values for k and 7 , and /q 0 t. The 
exceptions are those models in the degenerate families, with 
very large magnifications (WP94; WMS95) resulting from a 
profile which is flatter than isothermal. We note tentatively 
that Figs. |to@ show a weak correlation of confidence with 
magnification. These high magnification models will there¬ 
fore yield higher confidences due to the denser caustic net¬ 
work. The models of WP94 and WMS95, which have profiles 
significantly steeper than isothermal produce a total magni¬ 
fication within the range of our ensemble of macro-models, 
but with smaller k and larger 7 . 

The predicted Rba of the macro models (Tab. ^]) are 
marked on Figs. ^ and Most models predict an average 
flux ratio of between 0.8 and 1.1. From Fig. |(| this implies 
that, for a uniform optical source with Sopt = 0.025 rj 0 , the 
mid-IR source has Sir > lr/ 0 with > 90% confidence, and 
> 0.5 r/ 0 with > 95% confidence. The model of SWL98 yields 
an average flux ratio of ~ 1.0 — 1 . 1 , and hence an estimate of 
Sir > lrj 0 with > 97% confidence. A Gaussian source with 
a radius Sir = 2 a twice that of the assumed optical region 
is inconsistent with the results at > 99% for all but 2 models 
(K1 and I<5 for which we find limits of >98%). The mid-IR 
and optical flux may therefore be said to come from differ¬ 
ent emission scales with high confidence. The confidences on 
upper limits are far more systematically dependent on the 
theoretical flux ratio, ranging between 80% and 15% for the 
5.0 r/o source, and 95% and 30% for the 10.0ry o source. 
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Figure 7. Plots of source size limits as a function of mean flux ratio for different profiles. Panels (a)-(c): Lower limits ( Sir >0.25, 0.5 
and 1.0 rjo). Panels (d)-(f): Upper limits (Sir <1.0, 5. and 10.0?jo). The dots show limits for a circular disc source of radius ,S'/ r. The 
squares show results for an annular source with S ° ut = 2 and Sout = ^t-Sir. The optical source was uniform with Sopt = 0.025?7 o . 
Note that panel (d) = one minus panel (c). 


4.2 Different source profiles 

The above results assumed a Gaussian mid-IR source pro¬ 
file. However, the results should be dependent on the source 
profile assumed. If the mid-IR emission is due to dust, then a 
torus described by an annular profile may be a more appro¬ 
priate choice since there is a minimum radius inside which 
the dust cannot exist (e.g. AJB00). We have repeated the 


calculation of Eqn. |b| for all combinations of macro-models 
i — xi using circular top-hat profiles of radius Sir, and an¬ 
nular profiles with outer to inner radii ratios ^4^ = 2 and 

Sout = Sir (the factor of preserves flux for a uni¬ 
formly bright source). The left-hand panels of Fig. ^ show 
the resulting probabilities that Sir > 0.25 r] 0 , Sir > 0 . 5 r/ o 
and Sir > 1.0 r/ 0 as a function of the macro-model flux 
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ratio Rba■ Similarly, the right-hand 3 panels of Fig. ^ 
show the confidences that Sir < 1.0?; o , Sir < 5.0 r] 0 and 
Sir < 10.0r? o . The circular source limits are denoted by 
the dots, and the annular source limits by the squares. The 
consistency of the results for annular and circular profiles 
(particularly near (/is)/(/ m) = 1 ) suggests that the area of 
caustic structure covered by the projected source profile is 
the important quantity. This idea should extend to annuli 
of different inner to outer radii, and to the projected area of 
an inclined source. 


4.3 Other image combinations 

Figs. and ^ show a correlation between the confidence 
on mid-IR source size limits and Rba ■ This correlation was 
found for macro-models having k and 7 near the values 
found for images A and B by SWL98, but is also followed by 
the macro model parameters for image C (SWL98), which 
has values of k, and 7 that differ from image A by a factor of 
~ 2. This combined with the range of image magnifications 
over which the correlation holds leads us to assume that our 
ensemble of models can be used to describe the statistics 
for the other images as well. In particular, we use the en¬ 
semble to compute probabilities from Eqn. 6 for the other 5 
possible image combinations. Fig. ^ shows the resulting con¬ 
fidence that Sir > 1.0 r / 0 as a function of macro-model flux 
ratio Rij. The image flux ratios used were R°bai Racj Rda> 
Rbc, RRb and R°dci chosen so that R° PT ’° ba < Rif^ 
(note that only three of these ratios are independent). We 
used the optical fluxes from the 26th of September 1999. 
The predicted Rij of published macro-models (Tab. ^|) are 
marked on Fig. 

Strong limits on Sir are obtained from the flux ratios 
between images D and C, and B and C, as well as images 
B and A. In the former two cases, this is because all macro 
models predict R > 1. I 11 the latter case, images B and 
A yield the greatest agreement for Rij between different 
macro-models. These image ratios yield Sir > 1 ?; 0 with 

80% confidence for all macro-models. Rba and Rdc 
are independent ratios, and so their limits can be combined 
to produce a stronger limit of Sir > 1 ij 0 with 95% 
confidence for all macro-models. 


5 CONCLUSION 


Monito ring of the gravitationally microlensed Quasar 
Q2237+0305 has found significant microlensed optical vari¬ 
ability. Observations in July and September of 1999 showed 
(de-reddened) V-band flux ratios between the different sets 
of images that differed by large factors from the correspond¬ 
ing flux ratios in the mid-IR. In addition, the mid-IR flux 
ratios are similar to those measured in the radio, and to 
predictions of some lensing models, suggesting that mid-IR 
emission region is not subject to large microlensing varia¬ 
tion. The mid-IR emission region is therefore larger than 
the microlens Einstein Radius ( ij 0 ), and hence at least 2 or¬ 
ders of magnitude larger than the optical emission region, 
which is thought to be < 0.01»; o (Wambsganss, Paczyn- 
ski & Schneider 1990; Wyithe, Webster & Turner 2000). 
The colour difference between different images is due to the 
magnification/de-magnification of the optical emission. We 


have used microlensing models to calculate, as a function of 
mid-IR source size, the probability of obtaining the observed 
mid-IR flux ratios given the observed optical flux ratios. 

An alternative approach might be to note that the V- 
band flux in some images varied between the two mid-IR 
observations, while the mid-IR flux remained steady. In this 
case, the optical ratios need not be de-reddened, and as a 
result carry smaller errors. However, the variation in V-band 
is not significantly larger than the observational uncertainty 
in the mid-IR. We have therefore restricted our attention to 
flux ratio variation between different images. 

The flux ratio statistics for large sources are sensitive 
to the macro-model assumed. We have therefore computed 
source size limits for an ensemble of models. The proba¬ 
bility has the interesting property of being primarily sen¬ 
sitive to the macro-model flux ratio. As a result we have 
parameterised the macro-models available from the litera¬ 
ture by their predicted flux-ratio. The strongest limits are 
found from the flux ratios between images B and A. This is 
due in part to the differences between the relative V-band 
and mid-IR fluxes, as well as the more consistent macro¬ 
model predictions for these images. We find that the mid-IR 
source size Sir > lr) 0 with > 90% confidence, and > 0.5 r) 0 
with > 95% confidence. The IR-emission scale is larger than 
the optical emission with a confidence > 99%. 

The limit on the infrared source size derived here for 
the Einstein Cross may be converted to a limit on the 
brightness temperature. Assuming microlensing by stars, 
Sir <; r/ 0 ~ 10 17 cm, and the brightness temperature at 
10 /im (rest frame 3.7 fi m) is about Tb 5-j 7900K for a lu¬ 
minosity distance of 10 28 cm, magnification of 15, and flux 
of 20 mjy. This upper limit on the source brightness rules 
out non-thermal emission mechanisms, such as synchrotron, 
which typically have brightness temperatures of 10 8 ~ lo K. 
As argued in AJB00, the spectrum indicates thermal emis¬ 
sion at ~ 2000K, which is comparable to the sublimation 
temperature of dust. The flux from the QSO is sufficient to 
heat the dust at the sublimation radius of ~ 1 pc. 

One can turn this argument around: assuming that the 
IR emission is due to dust at the sublimation radius, then 
one can estimate the variance of the flux of each image for 
the given source size and macro lens model. We can then 
compute the y 2 of each macrolensing model using this vari¬ 
ance. Analytic estimates of the variance for large source exist 
for models with 7 = 0 ( Refsdal fe Stabell 1995) and some 
nume rical estimates exist for non-zero 7 (Refsdal & Stabell 


1997). The estimates with zero 7 indicate that the variance 


should be of order 10 % or less, unfortunately about the same 
size as the observational error bars for the infrared observa¬ 
tions. Since models for a large source size require much larger 
ray-tracing simulations, with a lower spatial resolution than 
we have carried out, we leave this computation for future 
work. Future observations with higher signal-to-noise, or a 
measurement of the variance with long-term IR monitoring 
will allow a better estimate of which macro-lens model is 
correct. 
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